A deep learning approach to 3D segmentation of brain vasculature

The segmentation of blood-vessels is an important preprocessing step for the quantitative analysis of brain vasculature. We approach the segmentation task for two-photon brain angiograms using a fully convolutional 3D deep neural network.Published versio

Cerebral tissue pO2 response to stimulation is preserved with age in awake mice

Published in final edited form as: Neurosci Lett. 2019 April 23; 699: 160–166. doi:10.1016/j.neulet.2019.02.007.Compromised oxygen supply to cerebral tissue could be an important mechanism contributing to age-related cognition decline. We recently showed in awake mice that resting cerebral tissue pO2 decreases with age, a phenomenon that manifests mainly after middle-age. To extend these findings, here we aimed to study how tissue pO2 response to neuronal stimulation is affected by aging. We used two-photon phosphorescence lifetime microscopy to directly measure the brain tissue pO2 response to whisker stimulation in healthy awake young, middle-aged and old mice. We show that despite a decrease in baseline tissue pO2, the amplitude of the tissue pO2 response to stimulation is well preserved with age. However, the response dynamics are altered towards a slower response with reduced post-stimulus undershoot in older ages, possibly due to stiffer vessel wall among other factors. An estimation of the net oxygen consumption rate using a modified Krogh model suggests that the O2 overshoot during stimulation may be necessary to secure a higher capillary O2 delivery to the tissue proportional to increased CMRO2 to maintain the capillary tissue pO2. It was observed that the coupling between the CMRO2 and capillary O2 delivery is preserved with age.Accepted manuscrip

Improved physiological noise regression in fNIRS: a multimodal extension of the General Linear Model using temporally embedded Canonical Correlation Analysis

For the robust estimation of evoked brain activity from functional Near-Infrared Spectroscopy (fNIRS) signals, it is crucial to reduce nuisance signals from systemic physiology and motion. The current best practice incorporates short-separation (SS) fNIRS measurements as regressors in a General Linear Model (GLM). However, several challenging signal characteristics such as non-instantaneous and non-constant coupling are not yet addressed by this approach and additional auxiliary signals are not optimally exploited. We have recently introduced a new methodological framework for the unsupervised multivariate analysis of fNIRS signals using Blind Source Separation (BSS) methods. Building onto the framework, in this manuscript we show how to incorporate the advantages of regularized temporally embedded Canonical Correlation Analysis (tCCA) into the supervised GLM. This approach allows flexible integration of any number of auxiliary modalities and signals. We provide guidance for the selection of optimal parameters and auxiliary signals for the proposed GLM extension. Its performance in the recovery of evoked HRFs is then evaluated using both simulated ground truth data and real experimental data and compared with the GLM with short-separation regression. Our results show that the GLM with tCCA significantly improves upon the current best practice, yielding significantly better results across all applied metrics: Correlation (HbO max. +45%), Root Mean Squared Error (HbO max. -55%), F-Score (HbO up to 3.25-fold) and p-value as well as power spectral density of the noise floor. The proposed method can be incorporated into the GLM in an easily applicable way that flexibly combines any available auxiliary signals into optimal nuisance regressors. This work has potential significance both for conventional neuroscientific fNIRS experiments as well as for emerging applications of fNIRS in everyday environments, medicine and BCI, where high Contrast to Noise Ratio is of importance for single trial analysis.Published versio

Development of a beam propagation method to simulate the point spread function degradation in scattering media

Scattering is one of the main issues that limit the imaging depth in deep tissue optical imaging. To characterize the role of scattering, we have developed a forward model based on the beam propagation method and established the link between the macroscopic optical properties of the media and the statistical parameters of the phase masks applied to the wavefront. Using this model, we have analyzed the degradation of the point-spread function of the illumination beam in the transition regime from ballistic to diffusive light transport. Our method provides a wave-optic simulation toolkit to analyze the effects of scattering on image quality degradation in scanning microscopy. Our open-source implementation is available at https://github.com/BUNPC/Beam-Propagation-Method.Accepted manuscrip

Contribution of speckle noise in near-infrared spectroscopy measurements

Near-infrared spectroscopy (NIRS) is widely used in biomedical optics with applications ranging from basic science, such as in functional neuroimaging, to clinical, as in pulse oximetry. Despite the relatively low absorption of tissue in the near-infrared, there is still a significant amount of optical attenuation produced by the highly scattering nature of tissue. Because of this, designers of NIRS systems have to balance source optical power and source–detector separation to maximize the signal-to-noise ratio (SNR). However, theoretical estimations of SNR neglect the effects of speckle. Speckle manifests as fluctuations of the optical power received at the detector. These fluctuations are caused by interference of the multiple random paths taken by photons in tissue. We present a model for the NIRS SNR that includes the effects of speckle. We performed experimental validations with a NIRS system to show that it agrees with our model. Additionally, we performed computer simulations based on the model to estimate the contribution of speckle noise for different collection areas and source–detector separations. We show that at short source–detector separation, speckle contributes most of the noise when using long coherence length sources. Considering this additional noise is especially important for hybrid applications that use NIRS and speckle contrast simultaneously, such as in diffuse correlation spectroscopy.R01 EB025145 - NIBIB NIH HHS; R24 NS104096 - NINDS NIH HHSPublished versio

Insight into the fundamental trade-offs of diffusion MRI from polarization-sensitive optical coherence tomography in ex vivo human brain

In the first study comparing high angular resolution diffusion MRI (dMRI) in the human brain to axonal orientation measurements from polarization-sensitive optical coherence tomography (PSOCT), we compare the accuracy of orientation estimates from various dMRI sampling schemes and reconstruction methods. We find that, if the reconstruction approach is chosen carefully, single-shell dMRI data can yield the same accuracy as multi-shell data, and only moderately lower accuracy than a full Cartesian-grid sampling scheme. Our results suggest that current dMRI reconstruction approaches do not benefit substantially from ultra-high b-values or from very large numbers of diffusion-encoding directions. We also show that accuracy remains stable across dMRI voxel sizes of 1 ​mm or smaller but degrades at 2 ​mm, particularly in areas of complex white-matter architecture. We also show that, as the spatial resolution is reduced, axonal configurations in a dMRI voxel can no longer be modeled as a small set of distinct axon populations, violating an assumption that is sometimes made by dMRI reconstruction techniques. Our findings have implications for in vivo studies and illustrate the value of PSOCT as a source of ground-truth measurements of white-matter organization that does not suffer from the distortions typical of histological techniques.Published versio

Colocalization of neurons in optical coherence microscopy and Nissl-stained histology in Brodmann’s area 32 and area 21

Published in final edited form as: Brain Struct Funct. 2019 January ; 224(1): 351–362. doi:10.1007/s00429-018-1777-z.Optical coherence tomography is an optical technique that uses backscattered light to highlight intrinsic structure, and when applied to brain tissue, it can resolve cortical layers and fiber bundles. Optical coherence microscopy (OCM) is higher resolution (i.e., 1.25 µm) and is capable of detecting neurons. In a previous report, we compared the correspondence of OCM acquired imaging of neurons with traditional Nissl stained histology in entorhinal cortex layer II. In the current method-oriented study, we aimed to determine the colocalization success rate between OCM and Nissl in other brain cortical areas with different laminar arrangements and cell packing density. We focused on two additional cortical areas: medial prefrontal, pre-genual Brodmann area (BA) 32 and lateral temporal BA 21. We present the data as colocalization matrices and as quantitative percentages. The overall average colocalization in OCM compared to Nissl was 67% for BA 32 (47% for Nissl colocalization) and 60% for BA 21 (52% for Nissl colocalization), but with a large variability across cases and layers. One source of variability and confounds could be ascribed to an obscuring effect from large and dense intracortical fiber bundles. Other technical challenges, including obstacles inherent to human brain tissue, are discussed. Despite limitations, OCM is a promising semi-high throughput tool for demonstrating detail at the neuronal level, and, with further development, has distinct potential for the automatic acquisition of large databases as are required for the human brain.Accepted manuscrip